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摘要 : 在 已 知 图 像 几何 扭曲 模型 的 基础 上 ， 对 2013 年 2 月 5 日 使 用 云南 天 文 台 2.4m 
望远镜 观测 近 地 小 行星 (99942 ) Apophis 资料 中 的 58 幅 CCD 图 像 进行 了 几何 扭曲 校正 ， 并 对 
校正 后 的 图 像 测 出 的 星 像 位 置 进行 了 几何 扭曲 残 差 分 析 。 结 果 发 现 ， 星 像 位 置 筷 、 了 两 个 方 
向 的 几何 扭曲 残 差 约 在 +0.02 pixel 内 。 为 了 探究 影响 几何 扭曲 残 差 大 小 的 因素 ， 设 计 了 仿 
真实 验 。 将 模拟 星 像 函数 经 过 设 定 的 扭曲 函数 模型 添加 几何 扭曲 ， 得 到 受 扭 曲 的 星 像 函 数 ， 
并 对 扭曲 前 后 的 星 像 进行 几何 扭曲 残 差 分 析 。 结 果 显 示 ， 星 像 的 半 峰 全 宽 (FWHM ) 越 大 ， 
几何 扭曲 模型 的 扭曲 变化 越 快 ， 几 何 扭曲 残 差 越 大 。 
关键 词 : 天 体 测 量 ; 几何 扭曲 ; 扭曲 校正 ; 仿真 
中 图 分 类 号 : P123.2*2 ”文献 标识 码 : A ”文章 编号 ; 1672-7673(2016)03-0382-07 


在 天 文 观测 中 ， 望 远 镜 光 学 成 像 系统 经 常 存在 几何 扭曲 ,导致 观测 得 到 的 天 体位 置 存在 几何 扭 
曲 误差 。 如 哈 勃 空间 望远镜 ( Habble Space Telescope, HST) 的 第 二 代 广 域 和 行星 照相 机 (Wide-Field 
Planetary Camera，WFPC2 ) 拍 摄 的 CCD 图 像 存 在 严重 的 几何 扭曲 。 其 非常 小 视 场 (WFPC2 视 场 中 的 宽 
场 仅 为 80 x 80 arcsec? ) 的 几何 扭曲 最 大 可 达 5 pixel。 这 种 由 光学 成 像 系统 导致 的 星 像 位 置 移动 ， 就 是 
几何 扭曲 ( Geometric Distortion, GD ) 。 文 [1] 通过 提高 线性 项 的 精度 提出 了 一 种 改正 几何 扭曲 的 方案 
求解 了 哈 勃 空间 望远镜 的 几何 扭曲 ， 使 得 哈 勃 空间 望远镜 拍摄 的 图 像 在 行星 相机 视 场 精确 到 0. 02 
pixel， 在 宽 场 视 场 精确 到 0. 01 pixel。 解 决 了 几何 扭曲 问题 后 ， 哈 动 空间 望远镜 才 真 正 发 挥 了 其 作为 
空间 设备 在 天 体 测量 上 的 潜能 。 在 地 面 望远镜 CCD 成 像 观测 中 ， 也 同样 存在 视 场 几 何 扭曲 的 影响 。 
文 [2-6] 对 地 基 天 文 望 远 镜 几何 扭曲 的 求解 做 了 系统 的 研究 ， 求解 了 云南 天 文人 台 2.4m 和 1 m 望远镜 
不 同 滤 光 片 的 几何 扭曲 模型 并 做 了 实际 应 用 。 本 文 基于 已 知 的 几何 扭曲 模型 改正 观测 的 CCD 图 像 ， 
并 对 其 中 产生 的 星 像 位 置 几何 扭曲 残 差 进行 了 仿真 研究 。 


1 观测 资料 和 几何 扭曲 模型 


1.1 观测 资料 

本 文 所 采用 的 图 像 是 2013 年 2 月 5 日 晚 使 用 云南 天 文 台 2.4 m 望远镜 配备 的 云南 暗 弱 天 体 光谱 
及 成 像 仪 (Yunnan Faint-Object Spectrog-raph and Camera, YFOSC) B 滤 光 片 观测 近 地 小 行 Apophis 所 得 
图 像 中 的 58 帧 。 因 为 云南 上 暗 弱 天 体 光 谱 及 成 像 仪 图 像 的 四 周 有 黑 框 ， 所 以 先 对 原始 图 像 进 行 适 当 的 
裁 甬 ， 裁 勇 后 的 图 像 大 小 为 1 900 x 1 900 像素 。 
1.2 几何 扭 遇 模 型 

文 [2-3] 提 出 了 一 种 简便 有 效 的 测定 望远镜 CCD 视 场 几何 扭曲 的 新 方法 ， 并 应 用 到 云南 天 文 台 
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中 星 像 的 定位 精度 有 了 明显 提高 ””。 本 文 使 用 


的 几何 扭曲 模型 正 是 用 这 种 方法 求解 得 到 的 。 Ns so. lobi IIT (1750,50) 1 (1850:50) | 
如 图 1， 将 大 小 为 1 900 x 1 900 像素 的 图 像 ”i '® 人 ? ? 
划分 为 19 x 19 个 大 小 为 100 x 100 RRITE + | : : HA] 
域 ， 黑 点 表示 子 区 域 的 中 心 像素 点 ( 格 点 ) 。 图 2 EOSO idir ERRI BOR 


为 求解 得 到 的 云南 天 文 台 2. 4 m 望远镜 B 滤 光 万 
的 几何 扭曲 矩阵 的 部 分 数值 。 和 矩阵 大 小 为 19 x 
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19， 其 中 每 一 个 数值 代表 对 应 格 点 位 置 的 几何 riit 
扭曲 量 。 对 于 该 几何 扭曲 模型 数值 对 应 的 精度 ， E (50.1750). | (18011750) (17530,1750) , (18530,1750) 
文 [2] 做 了 相关 研究 ， 结 果 显示 绝 大 部 分 几何 担 。 1000 o0 UT H ° 
E EA LBS boe 287] T Pope TH : (50,1850) ; (150,1850) (1750,1850) | (1850,1850) 
个 标准 差 的 分 布 是 相对 随机 的 ， 与 图 像 坐标 没有 io eco E E 
显著 的 相关 性 。 图 3 是 将 几何 扭曲 矩阵 放大 200 nca c ob e Ob A NCC R D CEDE unu) 
- 倍 后 用 矢量 方式 绘 出 。 从 图 3 可 以 看 到 ， 几 何 捏 图 1 格 点 
曲 在 图 像 的 上 面 角落 位 置 比较 大 。 mel Gak 
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图 2 几何 扭曲 模型 的 数值 矩阵 形式 图 3 几何 扭曲 模型 的 矢量 图 
Fig.2 GD pattern in matrix form Fig.3  Vectorgraph of GD matrix 


2 几何 扭曲 校正 


已 知 的 几何 扭曲 以 无 扭曲 位 置 为 参考 。 如 图 
4, Ñ P(xp, yp) 受 到 几何 扭曲 后 移动 到 已 (xmw， 
y,)， 记 为 PP'， 则 PP 点 的 几何 扭曲 为 dCP)= 
(dx(p), dy(p))， 如 (1) 式 。 

dx( P) 2xp — x, 
dy(P) =yp 一 xp 

申 己 知 的 几何 扭曲 矩阵， 图 像 中 任意 位 置 处 
的 几何 扭曲 可 由 离 该 位 置 最 近 的 4 个 格 点 的 几何 
扭曲 量 用 双 线 性 插值 的 方法 求 出 。 


(1) 


图 4 几何 扭曲 
Fig. 4 Geometric distortion 
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对 图 像 进 行 扭曲 校正 的 过 程 如 图 5。 黑 点 表 
示 一 幅 图 像 中 不 受 扭曲 的 整 像素 点 ， 受 到 几何 扭 
曲 后 偏 移 到 了 空心 点 的 位 置 。 如 整 像素 点 P(i， 
刀 受 到 扭曲 后 移动 到 已 (x, y) ， 那 么 用 (2) 式 求 
出 P' 的 位 置 。 在 扭曲 图 像 中 插值 求 出 已 位 置 的 
灰 度 值 ， 并 赋值 给 像素 P。 这 样 就 恢复 了 PP 像素 
点 无 扭曲 的 灰 度 值 。 对 图 像 中 的 每 一 个 像素 应 用 此 
过 程 ， 就 可 以 得 到 一 幅 扭曲 校正 后 的 图 像 。 
x =i + dx(p) 
y =j + dy(p) 


(2) 


Hp 
r 实验 结果 图 5 像素 扭曲 示意 图 
MER Fig.5 Distorted pixel 


对 观测 的 CCD 图 像 搜 星 得 到 存在 几何 扭曲 的 星 像 坐 标 ， 记 为 B(xu, ya)。 用 (3) 式 计算 其 扣除 
几何 扭曲 后 的 星 像 坐 标 C(x, y) 。 式 中 C 的 位 置 未 知 ， 其 扭曲 量 也 未 知 ， 所 以 无 法 直接 计算 C 的 位 
置 。 实 际 上 相 邻 位 置 的 几何 扭曲 几乎 相等 (实际 求解 的 几何 扭曲 模型 最 大 扭曲 量 约 为 +2 pixel) 。 所 以 
先 用 中 位 置 的 几何 扭曲 近似 C 位 置 的 几何 扭曲 ， 求 出 C 的 近似 位 置 ， 然 后 再 用 迭代 的 方法 计算 C 更 
精确 的 位 置 。 用 (4) 式 将 两 组 坐标 相 减 ， 得 到 星 像 的 几何 扭曲 。 图 6 中 灰色 点 为 星 像 位 置 的 几何 扭曲 
分 别 随 星 像 位 置 的 x、y 坐标 的 分 布 情况 。 

x 2x4 — dx(C) (3) | 
y 7 yq - dy CC) Ay 2 yq TY 

对 扭曲 校正 后 的 CCD 图 像 搜 星 得 到 的 星 像 中 心 ， 记 为 4(*xw,， yu). 与 C(x, y) 相 减 ， 得 到 扭曲 
校正 后 星 像 位 置 的 几何 扭曲 残 差 ， 即 (5) 式 。 图 6 中 黑色 点 为 星 像 位 置 的 几何 扭曲 残 差分 别 随 x、y 坐 
标的 分 布 情况 。 将 图 6 中 黑色 点 表示 的 几何 扭曲 残 差 量程 缩小 到 + 0. 05 pixel ， 得 到 图 7。 图 8 为 星 像 
的 几何 扭曲 残 差 放大 200 倍 后 的 矢量 分 布 图 。 


Ôx = Xa 一 区 


Ax =x —% 


(4) 


(3) 
ôy = ada =y 


* Xm * 立方 向 的 几何 扭曲 量 
* XIIL jk . ; * YALAR 


0 500 1 000 1 500 2 000 | 0 500 1 000 1 500 2 000 
X/pixel Y/pixel 


图 6 星 像 几何 捏 曲 和 几何 扭曲 残 差 随 X/Y 8 346 
Fig.6 The GD errors and residuals in X/Y 
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图 7 几何 扭曲 残 差 随 X/Y 的 变化 
Fig.7 The GD residuals in X/Y 


从 图 6 可 以 看 到 ， 图 像 做 扭曲 改正 前 ， 星 像 
的 几何 扭曲 最 大 可 达 1.7 pixel, 扭曲 改正 后 ， 星 M Dic LRL AL 
像 的 几何 扭曲 残 差 基本 分 布 在 0 附近 。 从 图 8 可 ou T ÉL ose tà ME ppt eth 
以 看 到 ， 星 像 的 几何 扭曲 残 差分 布 的 均值 为 
0. 008 pixel， 标 准 差 为 0. 009 pixel， 并 且 与 图 3 
相 比 明 显 消除 了 几何 扭曲 的 区 域 系 统 性 。 在 图 7 
中 将 残 差 量程 缩小 后 ， 可 以 看 到 几何 扭曲 残 差 大 
部 分 在 + 0. 02 pixel 内 ,但 发 现 几 何 扭 曲 残 差 的 
分 布 呈 现 一 定 的 系统 趋势 。 出 现 这 个 现象 ， 可 能 
的 一 个 原因 是 几何 扭曲 模型 在 图 像 边缘 处 的 变化 
比较 大 ， 造 成 图 像 改正 后 的 几何 扭曲 残 差 也 比较 
大 。 为 了 进一步 研究 影响 几何 扭曲 残 差 的 因素 ， 0 
本 文 做 了 星 像 几何 扭曲 的 仿真 实验 。 0 500 bud 1500 2000 

pixe 
4 仿真 实验 Es 几何 扭曲 残 差 失 量 图 
Fig.8 Vectorgraph of GD residuals 
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4.1 仿真 实验 设计 

如 图 9，(a) 为 用 一 维 连 续 高 斯 函数 模拟 的 星 像 f(x) ， 表 达 式 为 (6) 式 ,x 是 模拟 星 像 的 中 心 。 
(b) 为 用 正弦 函数 模拟 的 扭曲 函数 模型 g(x) ， 表 达 式 为 (7) 式 。 图 中 为 了 突出 星 像 扭曲 的 效果 ， 在 对 
扭曲 函数 绘图 时 缩小 了 扭曲 函数 的 周期 。 实 际 上 ， 一 个 星 像 范 围 内 的 扭曲 变化 非常 缓慢 。(e ) 为 对 模 
拟 星 像 函 数 添加 泊 松 噪声 ， 再 根据 扭曲 函数 添加 扭曲 后 得 到 的 连续 星 像 函 数 。 根 据 扭 曲 前 后 模拟 星 像 
函数 的 积分 面积 不 变 ， 由 (8) 式 可 以 推导 出 受 扭曲 后 的 连续 星 像 函 数 R(0) ， 如 (9) 式 。 


- (x-2)? 


2 
f(x)=B+He mm , (6) gla) = Asin T a), en 


x * g(x)-t _ Kx) 
f(x)dx 2 h(1)dt ' 3) aur XE 
根据 文 [7] 对 不 同 星 像 定 心算 法 的 比较 和 精度 分 析 表 明 ， 高 斯 拟 合法 是 一 种 精度 相对 较 高 的 定 心 
算法 ， 且 星 像 分 布 与 高 斯 函数 相似 程度 越 高 ， 则 测量 精度 越 高 。 实 际 求解 的 几何 扭曲 模型 最 大 扭曲 量 


(9) 
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心 。 具 体 地 ， 对 受 扭曲 后 的 连续 星 像 函 数 在 其 + 2. 50 像素 范围 内 采样 (o 是 高 斯 函数 的 标准 差 )。 设 
此 范围 内 覆盖 的 像素 区 间 为 (i, itn)。 即 对 (i, i+1), (itl, i+2)…(itn-1, itn) 区间 分 别 积分 ， 结 
果 得 到 一 组 离散 数据 点 ， 用 图 9(e) 中 的 点 表示 。 对 这 些 离散 数据 点 用 (10) 式 做 最 小 二 乘 拟 合 ， 求 出 
受 扭曲 后 的 星 像 中 心 x.'。 实 验 中 对 每 个 像素 区 间 积 分 时 ， 采 用 两 点 高 斯 数值 积分 ” 。 
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H(t) =B + H'e z7 . (10) 


(6)~(10) 式 中 , B 为 背景 值 ， 实 验 中 设 定 星 像 扭曲 前 后 背景 值 不 变 ; H 为 星 像 亮 度 值 ; 4 et K 
数 幅度 ， 实 验 中 参考 真实 资料 求解 得 到 的 扭曲 和 矩阵， 将 4 设 为 2 pixel; T 是 扭曲 函数 的 周期 ,7 越 小 ， 
HRR; H x o 为 对 离散 点 做 最 小 二 乘 拟 合 后 求 得 的 高 斯 函数 的 参数 。 
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图 9 (a) 模 拟 的 星 像 函数 ; (b) AH Ak; (c) 扭 曲 后 的 星 像 函数 


Fig.9 (a) Simulated star profile function; (b) GD function; (c) Distorted star profile function 


4.2 仿真 实验 结 
设 图 像 的 宽度 是 1024 pixel， 在 [0, 1023] 区 间 内 随机 生成 500 个 一 维 模 拟 星 像 函 数 。 由 (11) 式 计 
算 几 何 扭曲 残 差 ,实验 结果 如 图 9。 


ôd =x. -|x +g(x.)]. (11) 


图 10 中 (a) 表 示 设 定 扭曲 函数 的 周期 T=512 pixel， 模 拟 星 像 的 全 峰 半 宽 分 别 为 13、9、4 pixel 
时 ， 几 何 扭曲 残 差 随 星 像 * 坐标 的 分 布 情况 。 表 1 是 它 的 数值 统计 。(b) 表 示 设 定 模拟 星 像 的 FWHM = 
9 pixel， 扭 曲 函 数 的 周期 分 别 为 256、512、1 024 pixel 时 ， 几 何 扭曲 残 差 随 * 坐标 的 分 布 情况 。 表 2 
是 它 的 数值 统计 。 

从 图 10(a) 和 表 1 可 以 看 到 ， 当 扭曲 函数 的 周期 相同 时 ， 星 像 的 全 峰 半 宽 越 大 ， 几 何 扭 曲 残 差 
越 大 。 例 如 在 固定 周期 7=512 pixel 情况 下 ， 当 FWHM=4 pixel 时 ， 几 何 扭曲 残 差 最 大 为 0. 006 pixel, 
标准 差 为 0.001 86 pixel; 当 FWHM =9 pixel 时 ， 几 何 扭曲 残 差 最 大 增 至 0. 01 pixel， 标 准 差 增 大 到 
0. 002 76 pixel; ?4 FWHM = 13 pixel 时 ， 儿 何 扭曲 残 差 最 大 值 虽然 几乎 没有 增 大 , 但 标准 差 达 到 了 
0. 004 pixel, MBI 10(b) 和 表 2 可 以 看 到 ， 当 星 像 的 全 峰 半 宽 相 同时 ,扭曲 函数 的 周期 越 小 ， 几 何 扭 
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曲 残 差 越 大 。 例 如 在 固定 FWHM -9 pixel 情况 下 ， 当 7=1 024 pixel 时 ， 几 何 扭曲 残 差 最 大 为 0.007 5 
pixel; 当 T2256 pixel 时 ， 几 何 扭曲 残 差 最 大 达到 0. 015 8 pixel。 也 就 是 说 ， 星 像 的 全 峰 半 宽 越 大 ， 扭 
曲 函 数 的 周期 越 小 ， 几 何 扭曲 残 差 越 大 。 因 为 星 像 的 全 峰 半 宽 越 大 ， 扭 曲 郴 数 的 周期 越 小 ， 都 将 导致 
一 个 星 像 范 围 内 所 包含 的 扭曲 变化 越 大 。 从 图 10 还 可 以 看 到 ， 几 何 扭曲 残 差 的 变化 趋势 跟 扭曲 函数 
的 变化 一 致 。 图 7 中 几何 扭曲 残 差 的 分 布 在 图 像 的 四 周 比 较 大 ， 中 间 部 分 则 较 小 。 对 应 图 3 的 几何 扭 
曲 模型 ， 几 何 扭曲 在 图 像 的 角落 处 变化 较 快 ， 相 当 于 仿真 实验 中 的 扭曲 函数 的 周期 了 较 小 的 情况 ; 
几何 扭曲 在 图 像 的 中 间 部 分 变化 较 小 ， 相 当 于 仿真 实验 中 扭曲 函数 的 周期 了 较 大 的 情况 。 从 而 导致 
了 图 7 中 分 布 在 图 像 四 周 的 星 像 几何 扭曲 残 差 较 大 ， 分布 在 图 像 中 间 的 星 像 几何 扭曲 残 差 较 小 。 
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图 10 (a) 相 同 周期 下 不 同 全 峰 半 宽 的 扭曲 残 差分 布 ; (b) 相同 全 峰 半 宽 下 不 同 周期 的 扭曲 残 差 分 布 
Fig. 10 (a) The GD residuals in the same period but with different FWHM ; 
(b) The GD residuals with the same FWHM but in different periods 
表 1 相同 周期 不 同 全 峰 半 宽 下 扭曲 残 差 数值 统计 表 2 相同 全 峰 半 宽 不 同 周期 下 扭曲 残 差 数 值 统计 
Table 1 The GD residuals statistics in the same Table 2 The GD residuals statistics in the same 
period but with different FWHM FWHM but in different periods 
FWHM(T=512) Max SD T( FWHM =9) Max SD 
/ pixel / pixel / pixel / pixel / pixel / pixel 
4 0. 005 9 0. 001 86 256 0. 015 8 0. 005 64 
9 0. 010 1 0. 002 76 512 0. 008 7 0. 002 72 
13 0. 0103 0. 004 00 1 024 0. 007 5 0. 002 62 
+ » 
5 结 论 


在 已 经 求解 出 几何 扭曲 模型 的 情况 下 ， 对 CCD 图 像 进 行 几何 扭曲 校正 。 校 正 后 ， 星 像 位置 在 X、 
了 两 个 方向 的 几何 扭曲 残 差 约 在 + 0. 02 pixel 内 。 对 星 像 位 置 的 几何 扭曲 残 差 进行 仿真 研究 ， 分 别 测试 
了 星 像 在 相同 的 扭曲 孔 数 周期 下 不 同 全 峰 半 帘 ， 以 及 星 像 在 相同 的 全 峰 半 冤 下 不 同 扭曲 函数 周期 的 几 
何 扭 曲 残 差 的 大 小 和 变化 规律 。 结 采 发 现 星 像 的 全 峰 半 宽 越 大 ， 扭曲 函数 周期 越 小 即 扭曲 模型 变化 越 
快 ， 星 像 位 置 的 几何 扭曲 残 差 越 大 。 这 个 结论 解释 了 CCD 图 像 经 几何 扭曲 校正 后 分 布 在 图 像 四 周 的 
星 像 其 几何 扭曲 残 差 较 大 ， 而 分 布 在 图 像 中 间 的 星 像 其 几何 扭曲 残 差 较 小 的 现象 。 
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Residuals and Analysis of Geometric Distortion Correction of CCD 


Images Based on Geometric Distortion Pattern 
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(1. Department of Computer Science, Jinan University, Guangzhou 510632, China, Email; 8975332369 qq.com; 


2. Sino-French Joint Laboratory for Astrometry, Dynamics and Space Science, Jinan University, Guangzhou 510632, China) 


Abstract: Fifty eight CCD frames taken by the 2. 4m telescope of Yunnan Observatories on Feb 5, 2013 


were corrected for their Geometric Distortions ( GDs) with the solved GD pattern by the approach proposed in 


Reference |2]. It is found that the GD residuals of the pixel position of stars still exist about + 0. 02 pixels in 


both X and Y directions while they are bigger at some corners and smaller in the middle part of each CCD 


frame. A simulation experiment was then carried out to analyze the cause that affects the GD residuals. Our 
results show that the GD residuals grow bigger when the FWHM ( Full Width at Half Maximum) of a CCD 


stellar image becomes bigger or the GD pattern changes fast. This experiment can explain the variation of the 


GD residuals of stars after correcting the GDs. 
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